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Abstract 

We review some recent results on finite dimensional spin glasses by study- 
ing recent numerical simulations and their relationship with experiments. In 
particular we will show results obtained at zero and non zero temperature, 
focusing in the low temperature properties of the model, and contrast them 
with different pictures of the low temperature phase of spin glasses: Replica 
Symmetry Breaking, Droplet Model and Trivial-Not-Trivial Scenario. 



1 Introduction 

Spin glasses are still a problematic issue. The introduction of frustration and dis- 
order in a statistical model has posed a real challenge to both experimentalist and 
theoreticians. 

One can take, as an example, the "canonical spin glass": a metal (e.g. Copper) 
in which ferromagnetic impurities have been introduced (e.g. Manganese). This 
system can be studied in the RKKY framework and the result is an oscillating 
interaction which couples the magnetic moments of the material. This oscillatory 
behavior induced by the disorder (magnetic moments) also introduces frustration 
in the material [1]. 
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In this work, given the limitations of space, we have restricted ourselves to 
treating a few topics related to numerical simulations in finite dimensional spin 
glasses (only on Ising like models), focusing on the properties of the low temper- 
ature phase, yet we will treat them in detail. In the last years, a large amount of 
work on finite temperature numerical simulations but also zero temperature ones 
and experiments have been done. We will try to give a detailed description of 
some of these simulations and experiments, highlighting common observables in 
them and contrasting these results with some of the three main theoretical models: 
Droplet Model (DM) [2], TNT (Trivial overlap but Not Trivial link overlap) [3, 4] 
and Replica Symmetry Breaking (RSB) [5, 6, 7]. Related work can be found in [8]. 

Unfortunately, we have put aside in this work interesting studies on rheol- 
ogy [9, 10], ultrametricity [11, 12, 13], Heisenberg spin glasses [14, 15], suit- 
ability of the Edwards-Anderson model to describe real experiments [16], two di- 
mensional Ising spin glasses [17, 18], heterogeneity [19, 20, 21, 22], Sherrington- 
Kirkpatrick model [23, 24], sum rules [7], anisotropy [25], chaos [26, 27, 28, 29], 
eigenvalues analysis [30], three dimensional ferromagnetic spin glasses [31] and 
in field numerical simulations [32, 33, 34, 35]. The list of references given in this 
paragraph is not complete. 

Very good reviews and books have been written in the past years. We refer the 
reader to them [36, 37, 6, 1]. In addition numerical simulations have been reviewed 
in [39, 38], experiments in [40] and dynamics in [41, 42, 43, 44]. 

To put this work into context we will review (briefly) the three main theoretical 
approaches to spin glasses. 

The first one is the so-called Replica Symmetry Breaking. It is based in the 
standard procedure which has worked extremely well in Statistical Mechanics in 
the past decades (the paradigm is the ordered Ising model). Firstly, one must solve 
the model in the Mean Field approximation. This is equivalent to solve the infinite 
dimensional model exactly. This was done by G. Parisi in 1980 [5, 6]. His main 
results are that there exist a (countable) infinite number of (finite volume) pure 
states organized in an ultrametric fashion. The differences in extensive free energy 
among these pure states are of order one. In addition, the interface between two 
of these states is space filling (its surface scales as its volume, like an sponge) [7]. 
The Parisi solution also predicts a transition in magnetic field. 

Once we know which is the solution in infinite dimensions (where there are no 
fluctuations) we enable the system to fluctuate around the Mean Field solution (in 
this case, that of Parisi). The appropriate technique to handle this kind of problem 
is the Renormalization group (that can be implemented in the Field Theoretical 
approach) and the goals are computing the upper critical dimensions (above which 
Mean Field provides a good picture of the transition) and determining the criti- 
cal exponents (at fixed dimension or in the e-expansion) below the upper critical 
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dimension. Within this approach it is very difficult (since it is based mainly in 
perturbation theory) to estimate the lower critical dimension (the largest dimension 
below which there is not phase transition). The renormalization group program 
has been done (in part) by de Dominicis, Temesvari and Kondor [45]. We should 
remark that this approach does not change the low temperature properties. Hence 
in between the lower critical and infinite dimensions the qualitative description of 
the broken phase is still provided by the Parisi solution. 

Another compelling theory is the droplet model [2]. The rationale of this model 
is the Migdal-Kadanoff renormalization group. This technique is exact in one di- 
mension and is approximate in higher dimensions. The main results of the DM is 
that there are two pure states (only one, if we consider the global spin flip symme- 
try), and that the magnetic field destroys the phase transition. In addition we can 
mention that the typical excitations are compact domains of reversed (against the 
ground state) spins. The cost in energy of these excitations scales as a power of the 
typical size of the droplet, L 6 '. 

Recently has been proposed a third way which interpoles between the droplet 
model and RSB: the TNT proposal [3, 4]. In RSB 9 = 0, since we can create an 
excited state with 0(1) energy. In this third approach 9 = as in RSB but the link 
overlap is trivial (as in the droplet model: the probability distribution of the link 
overlap is delta peaked). In RSB the link overlap is believed to be proportional to 
the squared of the overlap (in infinite dimension the link overlap qi = q 2 , where 
q is the overlap): as far as the probability distribution of the overlap is not trivial 
then the probability distribution of the link overlap must not be delta peaked. We 
refer the reader to the text below for more details about the link overlap. 

The chapter is organized as follows. In the first section we examine the is- 
sue of the phase transition, giving strong numerical results which support a finite 
transition at non zero temperature. Next we will study the properties of the low 
temperature region (below the critical point). In this part we show numerical sim- 
ulations which highlight physical properties which can be described consistently 
assuming a RSB phase. Moreover we will describe the experimental computation 
of the dynamical correlation length and the possible interpretations of the different 
scalings proposed. In section 4 we will study the generalization of the fluctuation 
dissipation theorem out of equilibrium, starting with the analytical basis and con- 
tinuing with some numerical results which support the link statics-dynamics. This 
tool is very important because it can be implemented in experiments (we will show 
these). In section 5 we will show the memory/rejuvenations experiments. In the 
following section we will study zero temperature properties which probe the differ- 
ent theoretical pictures. Finally we will return to non zero temperature and describe 
recent numerical simulations computing the link overlap at finite temperature. 
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2 On the phase transition 



This part of the review is devoted to showing numerical evidences which favor 
strongly a phase transition in the three dimensional Ising spin glass at finite tem- 
perature. 

The existence of a phase transition in the three dimensional Ising spin glass 
has been attacked mainly using finite size scaling (FSS) methods [46]. In these 
methods one monitors which is the behavior of some (critical) observables of the 
system when one changes the size of it. We will describe in this section how to 
implementate the FSS to spin glasses and then how to define a good cumulant 
which signs clearly the transition point. 

The initial point is to introduce the Edwards- Anderson Hamiltonian [47] 

H = - ^2 Jij^j i (1) 

<i,3> 

where the sum is extended to all the pairs of nearest neighbors, <jj = ±1 are Ising 
variables and Jjj are random (quenched) variables. In general the Jij are drawn for 
a Gaussian distribution with zero mean and unit variance. One can also choose the 
random couplings from a bimodal distribution: = ±1 with equal probability. 

It is well known that observables in spin-glasses need to be defined in terms 
of real replicas, that is, for every disorder realization, one considers two thermally 
independent copies of the system {cij, Tj} [39]. Observables are most easily defined 
in terms of a spin-like field, the so-called overlap field (which is the order parameter 
in spin glasses): 

Qi = OiTi . (2) 
The total overlap is the lattice average of the % 

i 

while the (non-connected) spin-glass susceptibility is 1 



Xg = V(q 2 ) . (4) 

In Finite-Size Scaling studies, it is useful to have dimensionless quantities, that 
go to a constant value at the critical temperature. The standard example of this 
quantity is the Binder cumulant 



'As usual we use the brackets to denote the thermal average for a given choice of disorder, and 
the overline to mark the average over the disorder. 
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Figure 1 : g 4 and G cumulants as a function of the inverse temperature for the three 
dimensional ± J Ising spin glass [50]. 



Another example is the 52 cumulant [48], that measures the lack of self-averageness 
of the spin-glass susceptibility 



{q 2 } 2 _ {q 2 } 

92 = o • W 

(Q 2 ) 

In reference [51] a third cumulant was proposed which is a function of 54 and 52 

2 1 - 54 

These cumulants have been really useful to characterize phase transitions both 
in ordered systems (the Binder cumulant) and in disordered ones (52 in diluted 
Ising models). Nevertheless in Ising spin glasses they do not provide a clear sig- 
nature of the phase transition {i.e., a clear crossing between curves corresponding 
to different lattice sizes). In Figure 1 we show both cumulants as a function of the 
temperature for the three dimensional Ising model (with a binomial distribution for 
the couplings and helicoidal boundary conditions). 

Unfortunately, 54 and 52 require the evaluation of a four-point correlation func- 
tion, which is statistically a much noisier quantity than a two-point one. A more 
convenient observable is the correlation-length, which is defined only in terms of 
the two-point correlation function. Notice that its ratio with the lattice size is 
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again dimensionless [46]. We therefore are faced with the problem of defining 
a correlation-length on a finite lattice. This was done in Ref. [49]. The main steps 
of the constructions are the following. Let C(r) be the correlation function of the 
overlap field, 



C(r) 



1 

V 



(qm+r) 



(8) 



and C(k) its Fourier transform. Notice that (7(0) is the spin glass susceptibility. 
Then, inside the critical region on the paramagnetic side and in the thermodynam- 
ical limit, one has 



C(jfe) oc 



1 



k 2 + 2 
1 dC 



\k\\ «r\ 



C dk< 



(9) 
(10) 



On a finite lattice, the momentum is discretized, and one uses [49] a finite- 
differences approximation to eq. (10), 



e = 



i 



4 [sin 2 (fc*/2) + sin 2 (£&/2) + sin 2 (A^/2)] 



Xq 



C(fcm) 



(11) 



where \q was defined in eq. (4) and k m is the minimum wave-vector allowed for 
the boundary conditions used (e.g., k m = (2tv/L, 2tt/L 2 , 27r/L 3 ) for helicoidal 
boundary conditions). Of course, eq. (10) holds in the thermodynamic limit (L » 
£) of the paramagnetic phase. As we do not use connected correlation functions, £ 
has sense as a correlation length only for 2 T > T c . 

We can study the scaling behavior of the finite-lattice definition (11) on a crit- 
ical point, where the correlation function decays (in d dimensions) as r~( d ~ 2+v \ 
The behavior of the Fourier transform of the correlation function for large L in 
three dimensions is given by 



C(k) 



dr r 1 



_ r] sm(kr) 
kr 



(12) 



and one finds that Xq/C(k ra ) goes to a constant value, larger than unity, because 
ll^mll = 0(l/L). Furthermore, £/L tends to a universal constant at a critical point 
(like the Binder cumulant g±). Moreover, on a broken-symmetry phase, where the 
fluctuations of the order parameter are not critical, one has Xq = 0(L d ), while 

2 We use T c to denote the critical temperature obtained in numerical simulations and in theoretical 
computations and T s the one obtained in experiments. 
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Figure 2: £/L cumulant for the three dimensional ±J Ising spin glass (left part 
of the figure). The same plot for the two dimensional XY (ordered) model. Taken 
from reference [50]. 



C(k m ) = 0(1). Therefore the full description of the scaling behavior of £/L is as 
follows. Let be the correlation-length in the infinite lattice: in the paramagnetic 
phase, for L 3> Coo, one has £/L = 0(1/L). In the scaling region, where ^ > L, 
£/L = ( 1), while in a broken-symmetry phase on a lattice larger than the scale of 
the fluctuations, £/L = 0(L d / 2 ). Consequently, if one plots £/L for several lattice 
sizes as a function of temperature, the different graphs will cross at the critical one. 

We can see the (clear) crossing phenomena in Figure 2. Also shown in this 
figure (right part) the same observable for the two dimensional XY model (with 
no disorder). This double plot tells us that 1) there is a phase transition a finite 
temperature and 2) we should discard T c = and XY-like scenarios for the phase 
transition of the three dimensional Ising spin glass 3 . 

3 Some properties of the low temperature region 

In this section we will describe numerical simulations and experiments which try 
to discern which are the low temperature properties of the three dimensional Ising 
spin glass by working well below the transition point. We will start by discussing 

3 These numerical results have been obtained with the dedicated computer SUE [52], which has a 
performance of 0.2 ns/spin. 
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the properties at the upper critical dimension of the model (which is six) where 
the analytical predictions from RSB simplify. Then we will report results in three 
dimensions. Finally we will review some issues related to the behavior of the 
dynamical correlation length. 

3.1 d = 6 

First, we will check one of the RSB predictions. To do this, numerical work has 
been carried out just at the upper critical dimension. In this dimension there is 
no renormalization of the powers of propagators (i.e., the anomalous dimension 
vanishes) and only multiplicative factors occur 4 . 

If RSB holds, a l/p 4 propagator should be found by looking at the q = sector 
of the model in the broken phase (T < T c ) 5 . We remark that at six dimensions the 
equilibrium overlap-overlap correlation function constraint to q = was obtained 
by De Dominicis et al. [45] 

I T T <T e : u3) 

which corresponds to l/p 2 at T = T c (the usual critical propagator) and l/p 4 (the 
replicon mode) for T < T c . 

From this correlation function we can compute the associated (spin glass) sus- 
ceptibility 

X = Jd 6 x C(x) . (14) 

Since we are working on a finite lattice, the previous integral must be performed 
in a box of size L. If we want to observe the dynamical behavior of x, the upper 
limit in the integral should be changed to £ (t), the dynamical correlation length. At 
this point we can assume that ~ t l / z ^ T \ which defines an, in principle, effec- 
tive dynamical critical exponent, z(T). Furthermore, one can assume a functional 
dependence z(T) = 4T C /T, where T c is the critical temperature: with this temper- 
ature dependence we recover the value z at the critical temperature z(T c ) = 4 as 
predicted by Mean Field [6]. The result for x(t) is 

{t 1 / 2 if T = T 
!*ei if t<t c : (15) 

4 At the upper critical dimension logarithmic corrections appear. These have been studied numer- 
ically by Wang and Young [53], and subsequently computed analytically in reference [54]. For a 
discussion on the lower critical dimension see reference [55] 

5 This ergodic sector is very important. Out of equilibrium, the system remains in this sector. 
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This expression can be compactly written as 

x(t) ~ t h ^ 



(16) 



This formula should be valid if we remain all the time in the q = sector. Hence, 
the exponent h(T) is a discontinuous function of temperature: i.e., h(T~) = 1 
while h(T+) = 1/2. Moreover h(T), if the Ansatz for z(T) is right, should grow 
linear. 

This can be tested by performing an out of equilibrium numerical simulation in 
a large lattice. The run starts at random and suddenly the system is quenched below 
the critical temperature. At this point the growth of the non linear susceptibility is 
recorded. At the same time, one can check that the system (due to the large lattice 
simulated) develops no overlap (and so we are sure that we are simulating inside 
the q = sector of the theory). The strategy is to point out the discontinuity of 
the power of the q = propagator when we reach the critical temperature from 
below (the propagator changes from 1/p 4 to the standard and critical 1/p 2 prop- 
agator). So, one needs to redo the previous schedule but quenching to the critical 
temperature. 

In Figure 3 we plot the results and it is clear that the system behaves as RSB 
predicts: h(T) grows linear below the critical point and develops a discontinuity, 
just on the amount predicted by RSB, at the critical temperature. And so, it has been 
shown 1) the existence of the replicon mode at finite dimensions and 2) the growth 
of the correlation length can be described with the following law: £ (T, t) ~ 
with z(T) oc 1/T [56]. 

The next step is trying to see if this picture holds in lower dimensions, in par- 
ticular in the physical dimension three. 

3.2 d = 3 

In three dimensions, it is possible to handle this problem (replicon mode in addition 
to a given behavior of £(t, T)) by studying the decay with time and position of the 
overlap-overlap correlation function, 

C(X, t) = ^ (<Ti+xTi +x <7iTi)t ■ (17) 
i 

where a and r are two real replicas (which evolve with the same disorder) and 
the index i runs over all the points of the lattice. As usual we denote by (• • •) 
the average over the disorder and, in this context, ((• • -)) t is the average over the 
dynamical process (for a given realization of the disorder) at time t. In plain words, 
the two replicas (a and r) evolve with the same disorder but with different random 
numbers. 
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Figure 3: Susceptibility exponent as a function of the temperature in the six dimen- 
sional Ising spin glass [56]. Notice the linear region below the critical temperature 
and the discontinuity at the critical point. 
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Figure 4: Re-scaled correlation function, x a C(x,t), against the scaling variable, 
x/t 1 ^, for L = 24 and T = 0.35, 0.5 and 0.7. Taken from reference [75]. Notice 
the quality of the scaling. 
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Figure 5: We show four overlap-overlap correlations functions at T = 0.7. From 
top to bottom: 1) equilibrium correlation function, 2) equilibrium (but computed 
with a small cut-off on the overlap) correlation function 3) and 4) extrapolation 
to infinite time of two dynamical correlation functions computed using two dif- 
ferent annealing procedures. See the text for more details. Figure taken from 
reference [58]. 



In the q = sector (obtained simulating very large lattices, for large times, 
but by controlling that the overlap of the system is always very small) it has been 
obtained that the numerical data [59, 75] follow very well the following scaling 
law (it has been checked that this behavior also holds in four dimensions [57, 61]) 



We show in Figure 4 the scaling plot for three different temperatures and the 
fit using eq. (18). The scaling plot and the agreement with the fit is very good. We 
can cite that the a exponent does not show a clear temperature dependence in three 
dimension (a ~ 0.5) [59, 75] 6 , whereas in four dimension the situation is very 
different since the alpha exponent varies greatly with temperature [57, 61]. 

This scaling law provides us with the equilibrium form of the propagator, by 
taking the limit t — ► oo in eq. (18): 



6 In three dimensions [89] it has been found at zero temperature that a ~ 0.4, in good agreement 
with the value found at non zero temperature. 




(18) 



C eq (x) = Hm C(x,t)<x — 



t^oo X a 



(19) 
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where the proportionally constant is C eq (x = 1). Of course, the exponent is not as 
found in six dimensions due to the renomalization effects (see eq. (13)). 

A further test can be done in three dimensions in order to check the pure power 
law of the correlation function restricted to small overlaps, eq. (19). One can com- 
pare this behavior (obtained dynamically and using an extrapolation) with that ob- 
tained by computing at equilibrium the overlap-overlap correlation function by tak- 
ing only those measures with overlap q < 0.01. In Figure 5 we plot in the lower 
part of the figure two curves corresponding to the correlation function obtained in a 
dynamical process taking the extrapolation to infinite time. The upper curve is the 
equilibrium correlation function (computed without imposing cut-off) and finally 
the last curve is the equilibrium one computed using a small cutoff (g max = 0.01). 
The agreement between the lower three curves is really good. This plot provides an 
additional evidence to the existence of a replicon mode in three dimensions [58]. 
In the droplet model C cq (x) — ► g| A as x — > oo in contrast with the numerical 
results which support eq. (19). 

3.3 Dynamical correlation length 

We have seen that in six dimensions the correlation length can be fitted £(T, t) as 
t 1 / z ( T ) with z(T) = z c T c /T, where z c is the dynamical critical exponent at the 
critical point (T c ). In particular, it was found in three dimensions that [59, 75, 60, 
38] 

at,T)o,t°- 153 ^ T ^ , (20) 

where we have assumed that T c = 0.95(3). In four dimensions [57] a similar 
behavior was found 

Z(t,T)<xt - 19 W T / T ° , (21) 

where T c = 1.80(1). The behavior in four dimensions interpolates very well be- 
tween the three dimensional results and that obtained in six dimensions T) oc 

t 0.25T/T c _ 

This dependence of the dynamical correlation length with temperature and time 
has been checked experimentally. The basic idea of the experiment reported in ref- 
erence [62], was to introduce an external magnetic field and then operationally 
define the dynamical correlation length via the volume of the droplet which con- 
tributes to the Zeeman energy: -Ezccman oc N s xH 2 (where N s is the number of 
spins contributing to the Zeeman energy, H is the magnetic field and x is the mag- 
netic susceptibility). By effect of the magnetic field, the typical times of the dynam- 
ics are modified by a factor exp (-cN s xH 2 /T), where c is a numerical factor. By 
measuring this reduction factor one can extract the number of spins involved in the 
dynamics for a given waiting time and temperature and using that N s oc £(t w , T) 3 , 
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100 1000 



Figure 6: N s , number of spins participating in barrier quenching (and hopping) as 
a function of logi w at T = 0.78T g = 28 K for CuMn. The solid curve is the 
prediction for power law dynamics, while the dashed one comes from activated 
dynamics, see the text for more details. Taken from reference [62]. 



£(iw, T) can be computed. In this way they computed the correlation length and by 
performing the experiment at different temperatures. The following experimental 
dependence was found (see the solid line in Figure 6): 

/ f \0.169T/T g 

£(t w ,T) = 0.653 (^j = , (22) 

where 1/tq = 4.1 x 10 12 s _1 . The agreement with the result obtained in numerical 
simulations (see eq. (20)) is very good. Nevertheless a fit assuming activated dy- 
namics (droplet model) is also possible (see the dashed line in Figure 6), obtaining 

i 

0.21 

(23) 

However we see that the prefactor of the fit is really small (it would be natural 
for it to be 0(1)). Moreover the ip exponent is just at the lowest allowable value 
in the droplet model (ip = 0.2). However, numerical work suggests that 1/tp = 
1/0.7 [60, 38]. 

A plausibility argument for the linear dependence of the effective dynamical 
critical exponent z(T) with the inverse of the temperature was given by H. Rieger 
in reference [38]. Assuming an Arrhenius law and that free energy barrier for an 



f(t,r) = io- 



t (t w 

7^ log — 
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excitation of typical size L, scales as logL (i.e., tp = 0), we can obtain that 7 

r oc exp ( c-^^ I , (24) 



where c is a constant, which can be rewritten as r oc 

LV*CO withz(T) oc 1/T. 
However, a different argument based in the droplet picture and also accounting 
for the experimental data can be given [63]. Indeed, let us assume that the time 
needed to evolve a conformation on a scale of size l n is given by (this defines the 
ijj droplet exponent) 8 

/ ^pp \ 

t n = t(l n ) ~ r exp ( ^ J . (25) 

This behavior has been tested in Figure 6 and although the the fit is good the pa- 
rameters are not realistic enough (see above). Nonetheless it is possible to modify 
the previous formula in order to work in the neighborhood of the phase transition 

t n = t(l n ) ~ r Q Q exp (^f^j , (26) 

with Y(T) = To(T c — T) v ^. Near the phase transition this formula reduces to the 
usual (non activated) formula r ~ l n c . 

To test this generalization of the original droplet formula it is interesting to 
compute experimentally the following function (using the same procedure as in 
ref. [62]): 

_ fiog(W ; o)-W( t ..rn * 

In Figure 7, G(t w , T) is shown against T/T g for different waiting times, tempera- 
tures and three different spin glasses with different critical temperatures. The linear 
fit, supporting equation (26), is very good and the points extrapolate to near 1 when 
G approaches zero. 

In addition, Berthier and Bouchaud, in ref. [61], have tested this scenario via 
numerical simulations. In particular in four dimensions they have found that this 
droplet generalization works well. However the microscopic time they obtained in 
their fits shows (in three dimensions) a non monotonic dependence on the temper- 
ature, for which there is not physical explanation [61]. 

7 We introduce in this discussion the Boltzmann constant, fca, which has been set to one in the 
rest of the paper. 

8 By inverting in this formula l n in terms of t„ we obtain the activated dynamics prediction for 
the dynamical correlation length, see equation (23). 
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Figure 7: Plot of G(t w , T), defined in the text, against T/T g for different waiting 
times, temperatures and three spin glasses with different critical temperatures T g . 
The authors have used ip = 1.5, v = 1.3 and z c = 5 motivated by the experimental 
study of AgMn. If z c = 6 is assumed then the data extrapolate to T/T g = 1. The 
scaling is very good. From reference [63]. 



4 Off equilibrium fluctuation-dissipation relations 

One of the most important results of Statistical Physics at equilibrium is the so- 
called fluctuation-dissipation theorem. In this section we will review its theoretical 
basis and its generalization at early times in the dynamics. Moreover we will see 
how this generalization provides us with a useful tool to understand which are the 
properties of the low temperature phase at equilibrium. 

4.1 Theoretical basis 

The starting point is to perturb the original Hamiltonian, H, of a spin glass in a 
magnetic field as 

n' = n + J Ah(t)A(t) dt , (28) 

where 

Ti = - JijViVj + hJ2ai, (29) 

<i,j> i 

h being the magnetic field. We can define the following autocorrelation function 

C(h,t2) = (A{h)A(t 2 )) . (30) 
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Usually, A(t) = (Ji{t), and the associated response function 



SAh(t 2 ) 



(31) 

Ah=0 



The brackets (• • •) in eq. (30) and eq. (31) imply here a double average, one over 
the dynamical process and one over the disorder. 

In the dynamical framework, assuming time translational invariance, it is pos- 
sible to derive the fluctuation-dissipation theorem (FDT), that reads 

R(h,t 2 )=f39(ti-t2) dC t ,h) ■ ( 32 ) 

Ot 2 

where (5 = 1/T is the inverse temperature. 

The fluctuation-dissipation theorem holds in the equilibrium regime, but in the 
early times of the dynamics one expects a breakdown of its validity. Mean Field 
studies [66, 67, 68] suggest the following modification of the FDT (OFDR here- 
after): 

R(h,t 2 ) = pXiCfatiWh - t 2 ) dC( £ ,t2) • (33) 

ot 2 

where X defines the violation of fluctuation-dissipation. We can use the previous 
formula, eq. (33), to relate the observable quantities defined in eq. (30) and eq. (31). 
In the linear response regime, the magnetization can be written as (we report for 
completeness the formulas obtained for a Ising spin glass in a magnetic field, so 

m[h](t)^0) 



m[h + Ah](t) =m[h](t) 

1 



6h'(t') 



Ah{t') + 0{Ah 2 ) 



(34) 



h'(t)=h(t) 



and so, 



Am[h,Ah]{t) = f dt' R{t,t')Ah{t') + 0{Ah 2 ) , (35) 

where we have defined Am[h, Ah](t) = m[h + Ah](t) — m[h](t). Eq. (35) is 
just the linear-response theorem neglecting higher orders in Ah. By applying the 
OFDR we obtain the dependence of the magnetization with time in a generic time- 
dependent magnetic field (with a small strength), Ah(t), 

Am[h,Ah]{t) ~(3 f dt' X[C{t,t')} 8C ^P Ah{t') . (36) 

J —oo Ot 
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Next we let the system evolve with the unperturbed Hamiltonian until t = t w 
and then we turn on the perturbing magnetic field Ah (hence, the system feels a 
magnetic field h + Ah) 9 . Finally, with this choice of the magnetic field, we can 
write 

,dC(t,t') 



and 



Am[h,Ah](t) ~ Ah/3 ( dt' X[C(t,t')} X] } (37) 
Jt w of 

Am[h,Ah](t) ~ Ah/3 C duX[u], (38) 
JC(tU) 

where we have used the fact that we are working with Ising spins. In the equilib- 
rium regime (X = 1, as the fluctuation-dissipation theorem holds) we must obtain 

Am[h,Ah}(t) ~ Ah(3{\ -C(t,t w )) , (39) 

i.e., Am[h, Ah](t) T/ Ah is a linear function of C(t, t w ) with slope —1. 

In the limit t, t w — ► oo with C(t, t w ) = q, one has that X(C) — > x(q), where 
x(q) is given by 

x(q) = I" dq' P(q') , (40) 



where P(q) is the equilibrium probability distribution of the overlap with support 
[<7min, Qmax]- Obviously x(q) is equal to 1 for all q > q mSLX , and we recover FDT 
for C(t,t w ) > Qmax- This link between the dynamical function X{C) and the 
static one x(q) has been already verified for finite dimensional spin glasses [69]. 
The link has been analytically proved for systems with the property of stochastic 
stability [71]. 

We remark that we can use this formula to obtain g max as the point where the 
curve Am[h, Ah](t) against C(t, t w ) leaves the line with slope — (3 Ah. 
For further use, we define 

5(C) = Cdqx{q), (41) 
Jc 

or equivalently 



d?C 

C=q 

In the limit where I-tiwe can write eq. (38) as 



Am[Ah]{t)T 

— ~ b (C (t, t w )) . (43) 



9 In the first numerical application of this method, Franz and Rieger [70] chose another depen- 
dence of the magnetic field with time: h(t) = hoO(t w — t). 
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Looking at the relation between the correlation function and the integrated re- 
sponse function for large t w we can thus obtain q max , the maximum overlap with 
non-zero probability, as the point where the function S(C) becomes different from 
the function 1 — C. 

From the function S(C) we can get information on the overlap distribution 
function P(q), through eq. (42). Let us recall which is the prediction for the S(C) 
assuming the validity of each one of the competing theories described in the intro- 
duction. The droplet model predicts P(q) = 5(q — q) and consequently 10 

S{C) = { l OV C n ^V (44) 

v ' [ 1-C for C>g. 

On the other hand the RSB prediction for the overlap distribution^ ], P(q) = 
(1 - x M )S(q - ftnax) + x M S(q - ftnin) + p(q) (where the support of p(q) be- 
longs to the interval [g m i n , <Zmax]> <Zmin °c /i 4 / 3 and q max mainly depends on the 
temperature), implies that 



5(C) = 



imin j 



5(0) for C < 
s(C) for q min <C < q max , (45) 
1-C for C > q max , 



where s(C) is a quite smooth and monotonically decreasing function such that 

d 2 s{C) 



p(q) 



dC 2 



(46) 

C=q 



In Figure 8 we show three possible behaviors of the function S(C) (and for the 
closely related function P(q)). 

To finish this section we will recall an approximate scaling property of the 
probability distribution of the overlap that was introduced by Parisi and Toulouse 
(hereafter PaT) [72]. In particular in Mean field the PaT hypothesis implies 11 




10 In models with only one state, as the droplet model predicts for the Ising spin glass in a magnetic 
field, the equilibrium time is finite irrespective of the value of the volume of the system, hence, we 
can always thermalize any volume, and so the asymptotic behavior, for waiting times larger than the 
equilibration time, consists only of the straight line 1 — C. There is no horizontal part. 

"The goodness of this approximate Ansatz has been studied in reference [73] in the Mean Field 
approximation. They find that none of the Parisi-Toulouse scaling hypotheses about the q(x) behav- 
ior hold, but some of them are only violated at higher orders (taking as the parameter of the expansion 
the reduced temperature). 
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Figure 8: A possible model classification based on the function S(C). The big ar- 
rows represent delta functions. (A) corresponds to droplet model, (B) to one step of 
replica symmetry breaking and (C) to continuously broken replica symmetry (e.g., 
Parisi solution of an infinite dimensional Ising spin glass in absence of magnetic 
field). Taken from reference [82]. 



The result for C > q max is general (and true for finite dimension) and for 
<7min < C < q mSLX we make the following Ansatz: S(C) = AT(1 — C) B (in 
RSB A = 1 and B = 1/2). If we substitute this Ansatz in eq. (43) we obtain the 
following scaling equation 

H^ T -4> = /((!_ C )T-t>) , (48) 

where / is a scaling function and (f> = 1/(1 — B) (in Mean Field <fi = 2). In order 
to be consistent, the scaling function should be composed by a linear part (x) and 
by a power law part (Ax B ). 

In the rest of this section, we will discuss numerical simulations and experi- 
ments. 

4.2 Numerical Results 

In Figure 9 we show the numerical points obtained for two very large waiting times 
(in order to control that no dependence on t w is found) and the prediction from the 
statics: x(q). As a control we have computed the final point of the curve (the C = 
point) extrapolating at infinite time the magnetization using a power law fit. The 
agreement is very good. Notice that the asymptotic curve, which we can identify 
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Figure 9: Off equilibrium fluctuation-dissipation relation. We plot m(t,t w )T/h 
versus the spin-spin correlation function C(t, t w ) at T ~ 0.7T C . The lattice sim- 
ulated was L = 64 and we show two waiting times and two perturbing magnetic 
fields in order to control that linear response holds. The straight line with — 1 slope 
is the equilibrium prediction. All the points on this line are pseudo-equilibrium 
points. We have marked the extrapolation to infinite time of the susceptibility with 
the two leftmost points in the plot. Finally we have computed the x(q) function at 
equilibrium (from the numerical simulation, using parallel tempering [77, 78], of a 
L = 16 lattice). Taken from reference [69]. 



with the largest waiting time in the figure, is not compatible with the prediction for 
the droplet model (a horizontal part followed by the pseudo-equilibrium one) [69]. 

From this figure, we can compute the Edwards- Anderson order parameter (qEA = 
<?max) as the point at which the numerical points depart from the pseudo-equilibrium 
region (the straight line 1 — C). We can estimate qEA — 0.7. If the droplet model 
holds the order parameter should be qEA — 0.55 (in the DM, the asymptotic curve 
should be a horizontal straight line in the region [0, <tea]; tne fi na l point C = is 
provided by the infinite time extrapolation of the susceptibility, we can compute 
<?EA = 1 - mT//i| asyn ). 

We can test these possible values for qEA- To do this we recall equilibrium 
numerical simulation performed using parallel tempering [77, 78] in a wide range 
of lattice sizes: L = 4, 6, 8, 10 and 16. We plot in Figure 10 the equilibrium 
probability distribution of the overlap P(q). We can define qea(L) as the value of q 
in which P(q, L) shows a maximum. Furthermore we can analyze the dependence 
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Figure 10: Overlap probability distribution for L = 4,6, 8, 10 and 16 at T ~ 0.7T C . 
Taken from [58, 76]. 



of <Zea(£) with L. The simplest dependence is a power law: 

q EA (L)=q%> A +^- b , (49) 

where a and b are constants. In Figure 11, we show qea(L) versus L -L5 together 
with a linear fit [76]. Therefore, the data can be described with great accuracy 
assuming a power law with a non zero value of q EA ~ 0.7. Finally, the data does 
not support a power law fit with final value ~ 0.55. 12 Therefore we have obtained 
two compatible estimates of qEA at T = 0.7 using an off-equilibrium technique 
and an equilibrium one and both results agree in the statistical error. 

This technique can be implemented in experiments. This has be done in refer- 
ence [74] by studying the CdCri. 7^0.384 insulating spin glass with T g = 16. 2K. 
One measures the response and the autocorrelation between the spins. The first 
part of the work is not difficult, but the latter one has posed a challenge to the ex- 
perimentalists. We report in Figure 12 the plot of the violation of FDT In contrast 
with what happens in numerical simulations (where there is not a measurable de- 
pendence of the curves with the waiting time for the larger times simulated. See 
Figures 3 and 4 of reference [35] for a detailed study of the L and t w dependences), 

12 Incidentally, in the droplet model the probability of having an overlap different from the max- 
imum one (<7ea) goes to zero as L~ 6 . It is clear that data in Figure 10 rule out this possibility. In 
particular, P(0) — > const / 0. The same conclusion is reached if one works with window overlaps 
instead of the total ones [83]. 
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Figure 1 1 : Value of the overlap in which the probability distribution shows a max- 
imum as a function of the lattice size at T ~ 0.7T C . From [76]. 

in the experiment a strong dependence has been found for the reported curves with 
the waiting time, thus, an extrapolation to large (infinite) waiting time is manda- 
tory. This extrapolation is the dashed line shown in the figure. Notice also the 
dot-dashed line in Figure 12 which corresponds to the quasi-equilibrium regime. 
If one believes the extrapolation, the figure supports heavily the RSB scenario and 
discards that of the droplet model. 

Finally, we will end this section by showing a scaling analysis of the off- 
equilibrium fluctuation-dissipation relations. This has been done by using the PaT 
scaling which applies with great precision to the equilibrium probability distribu- 
tion, although it is not exact. In Figure 13 we report the scaling plot and it can 
be seen that is a really good scaling (different magnetic fields, in order to control 
linear response, waiting times, to check asymptoticity, and temperatures) [75]. 

Notice that the PaT scaling works for L and t w independent curves (see Figure 
9 and reference [35]). Two clear and distinctive regimes can be seen in that fig- 
ure. The first one correspond to the quasi-equilibrium regime: in that part of the 
figure the behavior is linear and thus it matches with the quasi-equilibrium regime 
AmT/h = 1 — C. The second one corresponds to the aging regime: that part 
of the plot can be parametrized with a power law with the B exponent introduced 
above. 13 

"Following reference [79] this kind of scaling is not enough to detect a RSB phase (they found 
in the two dimensional Ising model — with no phase transition at finite temperature — a PaT scaling 
for their OFDR). Nevertheless, in [79] the PaT scaling only works for points with the same waiting 
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Figure 12: Experimental determination of the function X which induces the vio- 
lation of fluctuation-dissipation. See the text for more details. Taken from refer- 
ence [74]. 
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Figure 13: PaT scaling for the three dimensional Ising spin glass (h = 0). The plot 
has been built with L and t w independent curves in order to check that we are in 
the asymptotic regime. Taken from reference [75]. 
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5 Memory and rejuvenation 



Maybe the most striking features of spin glasses were found in experiments where 
cycles in temperature were done. We are referring to the so-called memory and 
rejuvenation effects. 

In Figure 14 we show an experimental plot reported in Ref. [64]. In this exper- 
iment, a sudden quench from high temperature is done to a temperature (T = 12K) 
below the critical one (T g = 16.2 K). At this point the out of phase susceptibility 
is recorded. At a certain point of the experiment the temperature is lowered again 
(in this case to T = 10 K) and the out of phase susceptibility is recorded again. As 
can be seen in Figure 14 the out of phase susceptibility, measured at the new tem- 
perature (10 K), starts from a higher value than the susceptibility that the system 
had just before the quench from 12K to 10K. This is known as the rejuvenation of 
the system. When we cold a system it behaves as if it was younger than before, i.e., 
its out-of phase susceptibility is higher than the one the system had at the higher 
temperature just before the quench. In plain words, the system at the new, lower 
temperature is farther from equilibrium than in the last moments at the higher tem- 
perature. One can stay at the lower temperature for a while and then restore the 
temperature of the system to the original one {i.e., we heat the system from 10 K 
to 12 K). In Figure 14 we see that the system recovers the value of the out of phase 
susceptibility that it had just before it was cooled to 10 K. This phenomenon is 
known as memory effect. Notice in the inset of Figure 14 how, despite the strong 
relaxation produced at 10 K, the curves obtained in the higher temperature T = 12 
K in two separated time intervals are in smooth continuation. In reference [65] the 
reader can see good, recent and detailed experimental studies of rejuvenation and 
memory effects. 

Berthier and Bouchaud [61] have recently obtained rejuvenation and memory 
in the four dimensional Ising spin glass. We reproduce in Figure 15 their results. 
In three dimensions they have not seen these effects [61]. 

As a numerical approximation to the ac out of phase susceptibility Berthier and 
Bouchaud [61] proposed to use: 



where C(t, t') is the spin-spin dynamical correlation defined in eq. (30). 

time, instead, in the plot we have points computed with different waiting times. In effect, we remark 
again, the scaling reported in Figure 13 is t w -independent (at least in the numerical precision) which 
is a behavior completely different from the two dimensional spin glass (paramagnetic phase). For a 
paramagnetic phase and very long waiting time (i.e. all the points lie in the 1 — C straight line) the 
PaT scaling plot should consist in points over the linear part (quasi-equilibrium regime), and none in 
the power law part (aging regime). 




(50) 
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Figure 14: Rejuvenation and memory in a real Ising spin glass. Out of phase sus- 
ceptibility, x"( UJ ^a), of CdCr1.7Ino.3S4, with critical temperature T g = 16. 2K, 
during a cycle in temperature. The frequency, u, is 0.01 Hz and t a is the time 
elapsed since the quench. The figure inset shows that, despite the strong relaxation 
at 10 K, both part at 12 K are in perfect continuation of each other. From [64]. 
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Figure 15: Rejuvenation and memory in the numerical simulation of the four di- 
mensional Ising spin glass. Evolution of the ac correlation function following the 
schedule: T = 00 -> T x = 0.9 -> T 2 = 0.4 -> T\ -> T 2 T x . Notice that the 
critical temperature for this model is T c = 1.8. Taken from [61]. 
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Following Berthier and Bouchaud we can try to explain memory and rejuvena- 
tion effects in terms of the dynamical correlation length [61]. 

Let us first consider rejuvenation. The system at the higher temperature (T±) 
thermalizes its so-called fast modes (i.e., £(£,Ti) << L). When the system is 
frozen to a low temperature, these fast modes that are just equilibrated at the higher 
temperature are out of equilibrium in the new one, and so the system at the new 
lower temperature, is younger than before (£(£i, T2) << £(£1, T\), where £1 is the 
time the system elapses in T = T\). This mechanism does not rely on the concept 
of chaos in spin glasses [65] 14 . 

The memory effect can be understood as follows. We have said that rejuvena- 
tion involves the reorganization of small scales as compared to the lengths involved 
in the aging at T\, the higher temperature. When we heat the system from T2 to T\, 
these small scales "almost instantaneously" equilibrate at T\, and the aging restarts 
at T\ at the same point. More quantitatively: the time needed for the system to 
recover its age at T\ is given by ^memory, T\) ~ £(£2, T2), where £2 is the time 
elapsed in T^. If T\ — T2 is "large" then £ mcm0 ry << h- So, in this interpretation, 
memory is based in the existence of two, well separated, scales, while rejuvenation 
is based in the reorganization of small scales. 

6 Spin glass at zero temperature 

In the last years a large amount of numerical work has been devoted to numerical 
simulations at zero temperature. In particular has been studied the influence of 
perturbations in the ground state of the system. We will study in the next two 
subsections, two way to perturb the system and we will discuss the results in the 
light of the three scenarios (RSB, DM and TNT). 

6.1 Changing the boundary conditions 

We will review in this section numerical simulations performed at zero tempera- 
ture in which ground states of the system are computed with great accuracy (see 
reference [80]). 

14 Chaos, in this context, refers to the sensitivity of equilibrium states in the ordered phase to small 
changes in the couplings or in temperature. Temperature chaos postulates that typical equilibrium 
configuration at two different temperatures Ti and T2 respectively, are strongly correlated in a dis- 
tance lo which depends on AT = Ti — T2. For distances larger than lo the correlations between these 
two typical configuration go to zero. Scaling arguments provides lo ~ ATI 1 /™, where a = d s /2—6, 
d s is the fractal dimension of the interface (see next section) and 9 is the usual droplet exponent. The 
change of the equilibrium states as we change temperature explains rejuvenation. In Mean Field the 
effect of chaos in temperature is minimal [26]. In numerical simulations no clear chaos effects have 
been detected [27], but see [28]; for a droplet interpretation, see [29]. 
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One interesting observable is the link overlap (we will study in the next section 
its properties at finite temperature) defined by [81]: 



(51) 



where q(i) = a{i)T{i) is the overlap, where cr/r belong to a ground state that 
has been computed using periodic/antiperiodic boundary conditions respectively, 
and by /x we are denoting one of the d unitary vectors that can be defined in a 
d-dimensional hypercubic lattice (and so i + fj, is the point of the lattice neighbor 
of i in the direction provided by the vector /i). 

Neglecting the points at the boundary, we can define the interface between 
both ground state configurations as the region of space in which qi(i, /i) = —1. 
The probability to pick up such interface on a given random link, I, is given by 



where by qi we denote the disorder expectation of qi(i, y) averaged over all sites, 
i, and directions, fi. 

At this point one is faced with three possibilities (there are an additional fourth, 
but, we refer to the reader to reference [81] for a detailed explanation): 

1. The interface is confined to a region of width L z , with z < 1; inside of 
this region the interface could have overhangs. In this case p goes to zero 
following a pure power law L~ a , where a > 1 — z. 15 

2. The wandering exponent z is equal to 1 and the interface is a fractal object 
(not a multifractal) with fractal dimension, d s . Then p ~ L~ a with a = 



3. The exponent a is zero and the probability, p, goes to constant. Hence the 
interface is space filling. This last possibility is realized in the RSB scenario. 

Let us consider some examples. In a ferromagnet the ground state obtained 
with antiperiodic boundary conditions (a.b.c.) should be locally similar to that 
obtained with periodic boundary conditions (p.b.c), modulo an interface. As the 
interface, which is a flat surface, has no measure, the link overlap between these 
two ground states should be 1 in the large volume limit. In the droplet picture 
of spin-glasses the discussion is similar. The only difference with the previous 
example is that the interface is not necessarily flat rather it could be a corrugated 

15 Roughly, the volume of the interface is L d ~ x x L z . To obtain the probability we must divide 
the interface volume by the space volume, obtaining the relation between a and z. 



P = 




(52) 



d — d 
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surface (scenario 2). On the other hand, the RSB of spin-glasses is different: The 
ground state obtained under a.p.b. is expected to be locally similar to one of the 
low energy states of the spectrum of the p.b.c. case. Thus, the link overlap between 
these two ground states should tend in the large volume limit to a constant value 
different from 1 . 

By computing ground states, Marinari and Parisi reported the following values 
for the link overlap computed using four different methods [81]: qi = 0.755(15), 
0.80(6), 0.732(8) and 0.722(5) (these figures have been extrapolated to L -> oo). 
Putting all four results together, we can finally quote 



Note that qi is three standard deviations away from the droplet prediction qi = 1. 

We can try to recover this figure by performing numerical simulations at finite 
temperature and then try to extrapolate the data to zero temperature. This has been 
done in reference [84] by noting that since the overlap between the ground state 
computed with two different boundary conditions (since the change of boundary 
conditions can be regarded as a strong perturbation and the corresponding ground 
states are far away) is expected to be very small we can use the dynamical nu- 
merical simulations reported in the previous section (in which the overlap remains 
almost zero all the run) in order to obtain the value of the link overlap. It is easy to 
obtain qi since it is nothing but C eq (x = 1) (see eq. (19)). In Figure 16 we plot the 
values of C eq (l) obtained for different temperatures and we also mark the value 
obtained at zero temperature, see eq. (53). The consistency of both sets of data is 
very good. We will come back to this issue at the end of this section. 

6.2 Bulk perturbations 

Another way to perturb the system is to add a perturbation to the couplings of the 
model [4, 89]. The new Hamiltonian reads: 



where Nb = dL is the number of bonds, e is the strength of the perturbation and 
{er } is the ground state configuration computed with no perturbation (e.g., e = 0). 



qi = 0.79(7). 



(53) 




(54) 



The sum Y,<i,j> is extended over all pairs of nearest neighbors. 
We can define the following link overlap 



(a,0) \ ~> (10 „a n 

B <i,3> 



(55) 
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Figure 16: Values extrapolated to infinite time of the correlation function at dis- 
tance x = 1, Coo(l, T) = C eq (x = 1, T), versus T. We have also marked by two 
horizontal dotted lines the interval where the value computed in reference [81] lies, 
using T = ground state calculations. The consistency of the two results is clear. 



One overlap can be defined in the usual way 

i 

Hence, we can write eq.(54) as 

H' = H + eqf fi) , (57) 

where we have labeled the configuration {a} by the index a. 

One can show that the original ground state energy is shifted by an amount e 
(since g(°>°) = 1); the energy of any other state a (e.g., a low lying excitation of 
the spectrum) is shifted by an amount eq\ a,0 \ 

Let AE be the gap between the a state and the ground state one in absence of 
perturbation. If e > eqi + AE the new ground state should be the a one. 

One can compute the droplet prediction for the probability of this event. In the 
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DM the probability to have an excitation of energy E is given by 16 

P (E) = ^r/(^r) > (58) 

and so the probability to have an excitation with energy less than e(l — qi) is given 
by 



^< = I dE P{E) = g(-) , 



(59) 



where we have used that 1 — q% ~ L~( d ~ ds ^ and #(u) = J " (is /(s) 17 . Finally 
= 6' + d — d s . Now we can write the following scaling function for the overlap 
and the link overlap (which are given by the product of the probability of having 
a favorable droplet P < , times the contribution of this droplet to 1 — q and 1 — q[ 
respectively, see the previous footnote) 

1-9 « (60) 

We remark that by g and g7 we denote the average of the overlap and the link 
overlap over the disorder, respectively. For small e, assuming /(0) / 0, we obtain 
the following asymptotic formulas 

~ L d-d.+y ' (62) 
1 

~ L 2(d-d s ) + e> ■ ( 63 > 

Palassini and Young found for the three dimensional Gaussian spin glass with 
periodic boundary conditions [4] 

6' = 0.02(3) , d-d s = 0.42(2) , n = 0.42(3) , (64) 

and in four dimensions [4] 

0' = 0.03(5) , d-d s = 0.23(2) . (65) 

16 It has been assumed that the probability of an excitation on the ground sate scales as L 6 . In 
the droplet picture the 8 exponent is defined computing the scaling of the difference of free energies 
from boundary condition changes. It turns out that the droplet picture predicts 6 = 6'. 

In the droplet picture there are only two thermodynamic states related by spin-flip reversal sym- 
metry. The overlap between a droplet and the ground state is one minus a term which scales as the 
volume of the droplet (L d ) divided by the volume of the system (L d ). The link overlap in this cir- 
cumstance differs from 1 by a factor which is the volume of the interface of the droplet (L ds ) over 
the total volume (L d ). We are assuming that the excitation is one droplet with size proportional of 
that of the lattice size [4]. 
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This implies a non trivial probability distribution for the overlap but a trivial one 
for the link overlap. 

However the RSB scenario cannot be ruled out by the data because it is possible 
to fit 1 — q and 1 — qi to a constant plus scaling corrections, i.e., 

l-qi = a + ^ (66) 

with a = 0.28(3) (i.e., qi = 0.72(3)), b and c being positive constants. Indeed, 
these are scaling-corrections that are compatible with the RSB prediction: 9' = 
and d = d s . 

On the other hand, the droplets prediction is 6' = 9 ~ 0.2 (in three dimensions) 
and d — d s > 0. 

Krzalaka and Martin reached the same conclusions [3]. Next, Houdayer, Krza- 
kala and Martin [86] studied the topological properties of these excitations finding 
sponge-like conformations (this provides a geometrical picture [85] for the RSB 
scenario) which costs 0(1) in energy; finally they concluded that large finite size 
effects should be presented in order to explain the data with the RSB picture 18 . 
Related work by this group can be found in references [87]. 

Therefore they propose [3, 4], assuming the absence of strong scaling correc- 
tions, an intermediate or mixed scenario between droplet and RSB, the so called 
TNT picture (TNT for trivial (qi), non trival (q)). 

Yet, the controversy is not settled. In reference [89] Marinari and Parisi ana- 
lyzing the data assuming RSB obtained 

1 - qi (q = 0) = 0.245(15) (67) 

and by the study the correlation functions 1 — qi(q = 0) = 0.33(2). Here qi(q = 0) 
denotes that the link overlap has been computed using only configurations with 
mutual zero overlap. Moreover it has been found that qi (q) depends quadratically 
on q, as found in infinite dimension. 

In addition, reference [88] has obtained, by simulating the three dimensional 
Gaussian Ising spin glass but with free boundary conditions (in this reference the 
ground states were computed in an exact way) 

1-^ = 0.20(2), (68) 

where the superscript c denotes the average over those samples in which the unper- 
turbed and perturbed ground states are very different (i.e., the overlap is less than a 

18 Houdayer, Krzakala and Martin provide three values for the link overlap in their paper: q t = 
0.68, 9.72 and 0.75 depending on the number of parameters used in their fits. 
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given threshold value, g max = 0.4. Another cutoff q max = 0.2 gives essentially the 
same results). This value is in very good agreement with that obtained by Marinari 
and Parisi and cited above in this section (see eqs. (53) and (67)). However, the au- 
thors also found that if scaling corrections are allowed into the droplet prediction, 
then an equally good fit is found, and one obtains 

9' = 0.19(6) , d-d s = 0.44(3) , n = 0.63 . (69) 

Notice that 9' = 9 ~ 0.2 as the droplet picture predicts, but in contradiction with 
TNT. 19 



7 (More on) The link overlap (at finite temperature) 

The aim of this section is to study the properties of the link-overlap qi but at finite 
temperature [90, 91]. The goal here is to characterize the probability distribution 
of the overlap computing its variance. It has been found that 

var(gj) ~ L -w . (70) 

It is possible to compute the m exponent as a function of 9' and d s , the fractal 
dimension of an excitation. If we assume that this variance is dominated by the 
contribution of a single droplet of size L, then this event occurs with probability 
T/L e (assuming a constant density of states for these excitations, i.e., /(0) / 0, 
see eq. (58)). We have seen that 1 — q\ is proportional to L~( d ~ d °\ The same holds 
true for 5qi. Hence, the variance (which is the mean value of 5qf) is given by 

var(«) ~ ^L- 2 ^) , (71) 

and so \i\ = 9' + 2(d — d s ). The link overlap probability distribution obtained 
with numerical simulation of the three dimensional Ising spin glass with periodic 
boundary conditions is shown in Figure 17. In the droplet model this probability 
distribution should shrink to a Dirac delta (zero variance); instead, in RSB P(qi) 
should have a compact support (and so non zero variance). 

The extrapolation of the /// exponent to zero temperature gives a value: [i\ = 
0.76(3). By assuming 9' = this implies that d — d s = 0.38(2), near the value 
computed directly at T = using ground state computations: d — d s = 0.42(2) 
(see preceding section). However, a fit assuming RSB (i.e., fi' = 0) cannot be ruled 
out by the numerical data [90]. 

''incidentally, the relation between the link overlap and the overlap has been study in this 
reference [88] obtaing qi = 0.77(2) + 0.27(3)<? 2 , according with RSB predictions. Moreover 
qi(q = 0) =0.77(2). 
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Figure 17: P(<ft) for the three dimensional Gaussian spin glass at T = 0.2 for dif- 
ferent lattice sizes with periodic boundary conditions. Taken from reference [90]. 



In reference [91] free boundary conditions (f.b.c) were used. With this f.b.c 
they try to discern between a trivial behavior for the link overlap given by var(gj) ~ 
L _e with a suitable e exponent, and the RSB behavior given by vax(qi) = a+bL~ c . 
They found a finite value for a, for all the temperatures simulated, which implies 
d = d s and that a pure law behavior (L~ e ) is excluded by the data. In Figure 18 
the variance of the link overlap is plotted against the lattice sizes, in addition to 
the different fits used (to a + b/L c ). The same analysis on the four dimensional 
Ising spin glasses provides the same picture (d = d s and 9' = 0). However, let 
us end this section recalling that it has been argued [88] that, in principle, results 
obtained with f.b.c show larger finite-size corrections than results obtained with 
p.b.c. On the other hand, f.b.c do not pose any restriction on the position of the 
domain wall. Hence, it is not clear if the results reported for f.b.c represent the 
asymptotic behavior and what are the optimal boundary conditions for this kind of 
studies. 

8 Conclusions 

We have reviewed in detail (some) recent works on finite dimensional spin glasses 
obtained with numerical simulations and in some cases, we have done a direct 
comparison with experimental results. Moreover we have tried to describe the 
numerical data with the competing three theoretical pictures which try to describe 
the low temperature phase of finite dimensional spin glasses. 
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Figure 18: Log-log plot of the variance of the link overlap as a function of the 
size (with free boundary conditions) for different temperatures. Taken from refer- 
ence [91]. 



Of course, as said in the introduction, we have omitted important issues (we 
apologize), but we hope that this review will clearly show which are the difficulties 
faced and the open/closed problems regarding this interesting and active field of 
Statistical Mechanics and Condensed Matter. 
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